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A  faster  method  for  computing  optimal  three  dimensional  tra¬ 
jectories  that  maximize  the  landing  footprint  of  a  lifting  re-entry 
vehicle  has  been  developed.  The  method  utilizes  energy  approxima¬ 
tions  based  on  the  assumption  that  the  flight  path  angle  is  small 
and  the  flight  path  angular  rate  is  zero.  Thus,  the  vertical  com¬ 
ponent  of  lift  is  considered  equal  to  the  weight  minus  centrifugal 
relief,  and  the  equation  of  motion  are  reduced  in  order  from  six 
to  four.  Because  of  this  simplification,  the  classical  indirect 
method  of  the  calculus  of  variations  is  used  to  compute  families 
of  optimal  solutions.  Utilizing  data  corresponding  to  one  of  the 
space  shuttle  configurations,  computations  have  been  carried  out 
for  both  unconstrained  trajectories  and  for  solutions  that  have 
aerodynamic  heating  rates  and  lift  coefficients  limited  to  speci¬ 
fied  values.  During  the  investigation,  several  interesting  ana¬ 
lytical  finds  were  uncovered  that  could  be  used  as  a  basis  for  an 
onboard  guidance  scheme. 
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ABSTRACT 

A  faster  method  for  computing  optimal  three  dimensional  tra¬ 
jectories  that  maximize  the  landing  footprint  of  a  lifting  re-entry 
vehicle  has  been  developed.  The  method  utilizes  energy  approxima¬ 
tions  based  on  the  assumption  that  the  flight  path  angle  is  small 
and  the  flight  path  angular  rate  is  zero.  Thus,  the  vertical  com¬ 
ponent  of  lift  is  considered  equal  to  the  weight  minus  centrifugal 
relief,  and  the  equations  of  motion  are  reduced  in  order  from  six 
to  four.  Because  of  this  simplification,  the  classical  indirect 
method  of  the  calculus  of  variations  is  used  to  compute  families 
of  optimal  solutions.  Utilizing  data  corresponding  to  one  of  the 
space  shuttle  configurations,  computations  have  been  carried  out 
for  both  unconstrained  trajectories  and  for  solutions  that  have 
aerodynamic  heating  rates  and  lift  coefficients  limited  to  speci¬ 
fied  values.  During  the  investigation,  several  interesting  ana¬ 
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INTRODUCTION 


The  introduction  of  energy  approximation  methods  (Refs.  1-5) 
in  1944  has  taken  the  science  of  high-speed  aircraft  performance 
calculations  out  of  the  realm  of  quasi-steady  state  analysis  and 
into  the  framework  of  dynamics  and  the  calculus  of  variations. 

These  early  attempts  of  extreme  reduced-order  approximations  were 
applied  only  to  two  dimensional  minimum  time  to  climb  problems  and 
resulted  in  solutions  that  were  plagued  with  subarcs  consisting  of 
unrealistic  climbs  and  dives  and  large  instantaneous  flight  path 
changes.  Nevertheless,  these  solutions  did  resemble  actual  opti¬ 
mum  flight  paths  and  provided  pilots  with  procedures  for  improved 
aircraft  performance. 

During  the  past  twelve  years  the  use  of  optimization  tech¬ 
niques  developed  specifically  for  modem  digital  computers  (Refs.  6 
10)  has  gone  to  the  other  extreme,  that  of  generating  numerical  so 
lutions  to  complex  three  dimensional  trajectory  problems  involving 
specified  terminal  conditions  and  inequality  constraints  on  control 
and  state  variables.  Although  these  solutions  are  realistic  and 
provide  a  good  reference  of  maximum  possible  performance,  it  is 
difficult  to  obtain  an  insight  to  the  general  characteristics  of 
optimum  aircraft  maneuvers,  and  the  high  cost  of  computer  usage 
makes  it  prohibitive  to  carry  out  any  extensive  aircraft  parametric 
design  studies.  It  would  be  desirable  to  have  available  a  simpler 
and  more  efficient  method  that  could  overcome  the  preceding  diffi¬ 
culties  . 

One  promising  approach  that  has  been  revived  and  extended  is 
the  energy- state  or  energy-maneuverability  approximation  (Refs.  11- 
14)  .  The  essential  feature  of  the  energy  approximation  is  the  re¬ 
duction  in  order  of  the  differential  equations  of  motion,  and 
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therefore  a  corresponding  reduction  in  the  Euler  equations  and  the 
number  of  multiplier  initial  values.  Although  the  reduction  may 
not  be  sufficient  for  deriving  analytic  solutions  it  should  facili¬ 
tate  use  of  the  classical  indirect  method  of  the  calculus  of  varia¬ 
tions  for  efficiently  computing  many  families  of  optimal  solutions. 
This  approach  has  been  successfully  carried  out  recently  in  a 
series  of  studies  for  calculating  optimum  turning  maneuvers  of  air¬ 
craft  at  constant  altitude  (Refs.  15-18)  and  in  three  dimensions 
(Refs.  19-22). 

The  research  investigation  described  in  this  report  extends 
the  energy  type  of  approximation  to  the  computation  of  optimum 
three  dimensional  trajectories  that  maximize  the  landing  footprint 
of  a  lifting  re-entry  vehicle.  The  footprint  may  be  considered  as 
a  locus  of  points;  for  each  point  the  terminal  condition  is  speci¬ 
fied  and  the  crossrange  for  any  value  of  downrange  has  been  maxi¬ 
mized.  The  approximation  is  based  on  the  assumption  that  the 

flight  path  angle  y  is  small  (cos  7*1)  and  the  angular  rate 
• 

y  is  zero.  Thus,  the  vertical  component  of  lift  is  equal  to  the 
weight  minus  centrifugal  relief.  This  has  the  effect  of  elimi¬ 
nating  the  skipping  type  of  trajectory,  and  instead  results  in  an 
equilibrium  type  of  optimum  glide.  It  should  be  noted  that  it  is 
the  long  duration  skipping  motion  that  plays  havoc  with  the  conver¬ 
gence  of  iterative  optimization  techniques. 

In  addition  to  maximizing  the  landing  footprint,  the  computer 
program  also  satisfies  inequality  constraints  on  underbody  tempera¬ 
tures  and  maximum  lift  coefficient.  No  attempt  has  been  made  to 
optimize  trajectories  that  terminate  inside  the  footprint;  however, 
this  could  conceivably  be  achieved  by  decreasing  the  value  of  the 
maximum  allowable  underbody  temperature. 
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Because  of  the  reduced  order  of  the  system,  there  is  a  problem 
of  transitions  from  and  to  specified  end  conditions.  As  in  the 
earlier  minimum  time  to  climb  problems  (Ref.  4),  where  energy  is  a 
state  variable  and  altitude  (or  velocity)  is  employed  as  a  control 
variable,  a  discontinuity  in  the  control  variable  is  mathematically 
permissible  and  results  in  instantaneous  altitude/velocity  changes. 
What  is  needed  to  correct  these  unrealistic  maneuvers  during  re-entry 
is  a  method  rf or  getting  on  and  off  the  optimal  equilibrium  glide  path 
in  an  optimal  manner.  Since  the  reduced-order  modeling  eliminates 
the  higher  frequency  motion,  these  fast  transients  have  to  be  re¬ 
introduced  and  superimposed  on  the  slowly  varying  motion  in  order 
to  satisfy  the  boundary  conditions.  These  transitions  are  suitable 
for  treatment  by  singular  perturbation  theory  and  the  "boundary 
layer"  correction  given  as  solutions  of  differential  equations  also 
of  reduced  order  (Refs.  13  and  22-24).  Some  preliminary  work  along 
these  lines  has  been  carried  out,  but  is  not  included  in  this  re- 


MATHEMATICAL  MODEL 


A  computer  program  has  been  developed  to  optimize  the  three 
dimensional  turning  flight  of  a  lifting  re-entry  vehicle,  taking 
into  account  constraints  placed  on  underbody  temperatures  and  maxi¬ 
mum  lift  coefficient.  A  spherical  nonrotating  earth  with  an  inverse 
square  gravitational  field  and  an  exponential  atmosphere  is  assumed. 
Rigid  body  dynamics  are  neglected,  and  throughout  the  flight  the 
vehicle  executes  a  coordinated  turn  maneuver  of  zero  side  force. 
Angle  of  attack  and  bank  angle  are  the  two  control  variables  for 
the  original  model;  however,  for  the  reduced-order  model,  altitude 
replaces  angle  of  attack  as  a  control  variable.  The  aerodynamic 
data  are  independent  of  Mach  number. 


Equations  of  Motion 

The  flight  path  equations  before  reduced-order  approximations 
are  made  have  the  form 
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where  V  is  velocity,  7  flight  path  angle,  x  heading  angle. 


r  geocentric  radius,  cp  latitude,  A  longitude  (see  Figs.  1  and 

2  2 

2)  ,  and  the  gravity  g  -  gQro/r  .  These  equations  are  derived  in 
Appendix  A. 

The  specific  energy  E  replaces  velocity  through  the  defining 
relationship: 


The  choice  of  specific  energy  (foot-pounds  per  pound)  as  a  state 
variable  stems  from  the  intuitive  reasoning  that  E  should  be  more 
slowly  varying  than  either  V  or  h.  This  would  be  particularly 
true  if  reduced-order  approximations  were  not  used,  resulting  in  a 
skipping  motion  that  can  be  thought  of  as  an  interchange  of  kinetic 
and  potential  energy  about  some  equilibrium  energy  level.  There  is 
another  reason  for  the  preference  of  E  over  V.  Because  E  de¬ 
creases  monotonically  it  is  possible  to  reduce  the  order  of  the 
system  by  one  by  eliminating  time  and  employing  E  as  the  indepen¬ 
dent  variable.  As  shown  in  Appendix  B,  the  order  of  the  system  is 
further  reduced  from  four  to  three,  thereby  resulting  in  an  analyti¬ 
cal  solution  for  the  unconstrained  case. 


Substitution  of  Eq.  (4)  into  (1)  eliminates  the  sin  7  terms 
from  the  equations  of  motion,  and 


V 


(8) 


This  is  significant,  since  the  small  angle  approximation  to  be  made 
for  7  applies  only  to  the  cos  7  terms,  i.e.,  cos  7  «  1.  The 
sin  7  terms  do  not  have  to  be  approximated  by  sin  7  w  0  in  order 
to  derive  exaoi  expressions  for  E  and  E.  Differentiating  Eq.  (7) 
and  combining  with  Eq.  (8)  results  in 
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and  V  has  been  reduced  to  an  ordinary 


where  the  weight  W  ■  mgQ, 
variable  defined  by 


v2 


2goE  + 


2 

2g  r 
o  o 


(9) 


Based  on  the  assumption  that  cos  y  *  1  and  the  angular  rate 
7*0,  Eq.  (2)  is  converted  to  an  algebraic  expression 
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Vo 
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where  the  altitude,  h  *  r  -  rQ,  is  now  employed  as  a  control  vari¬ 
able  in  addition  to  the  angle  of  attack  a  and  bank  angle  p..  Ihe 
system  of  equations  is  reduced  to  fourth  order: 


E 


DV 

W 


•  V  Sin  ^  V 

X  *  - — - -  sin  v  tan 

•  V 

cp  =  -  —  COS  Y 
r 

•  V  sin  y 
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(U) 
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(13) 

(14) 


Because  of  Eq.  (10)  the  three  control  variables  are  not  inde¬ 
pendent  of  each  other.  Instead,  h  and  p.  are  to  be  considered 
as  the  two  primary  control  variables  for  the  reduced-order  system 
and  a  is  to  be  determined  from  Eq.  (10)  and 


p  - 

-Ph 

po  6 

(15) 

CL- 

L 

(16) 

}p\rs 

a  « 

.145  +  .737 

tan (1.33  CL) 

(17) 

where  the  numerical  expression  for  a  versus  C^,  with  lim¬ 

ited  to  0.535,  corresponds  to  one  of  the  early  space  shuttle 
orbiter  configurations.  Other  equations  and  data  are: 


(18) 


D 


cDlpv2s 


(19) 


u 

LI 
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where  CD 

o 


1.69,  W  -  211,170  lbs. 


PQ  *  .0034  slug/ft"*,  p  ■  .0000421  per  ft,  rQ  ■  20,908,800  ft, 
gQ  =  32.174  ft/sec2  and  S  =  6100  ft2. 


Boundary  Condi t ions 

At  the  start  of  re-entry  Eq  =  -  11,717,000  ft,  yq  *  tt/2  radians 
(heading  east),  <pq  =  0  (on  the  equator)  and  Aq  *  0.  The  initial 
condition  for  E  corresponds  to  a  velocity  of  24,000  ft/sec  at 
an  altitude  of  40  nautical  miles.  However,  since  altitude  is  a 
control  variable,  to  be  determined  optimally,  only  Eq  remains 
fixed  for  each  member  of  the  footprint  family,  whereas  the  numeri¬ 
cal  values  of  altitude  and  velocity  at  t  ■  0  will  vary  while 
satisfying  Eq.  (7). 
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The  fln.il  conditions  of  the  turning  flight  re-entry  are 
Ej.  =  -  20,800,000  ft,  Yp  is  open  and  rip  and  Ap  are  arrived  at 
optimally  such  that  the  area  of  the  landing  footprint  has  been  maxi- 
mi  zed.  Hie  value  of  Ep  corresponds  to  a  velocity  of  2000  ft/sec 
at  an  altitude  of  50,000  ft. 

tnequa I i iy  Cons t r a ints 

The  underbody  temperature  constraint  is  taken  from  Ref.  25, 
and  for  this  investigation  is  considered  to  be  typical  for  those 
cases  where  the  boundary  layer  is  assumed  to  be  in  thermal  and 
chemical  equilibrium.  The  constraint  is  given  as  boundaries  in 
al titude- velocity  space  (Fig.  3)  such  that  all  underbody  tempera¬ 
tures  aft  of  the  20-foot  station  are  limited  to  1800  F.  At 
higher  altitudes  the  maximum  temperature  would  be  less  than  1800  F 
and  at  lower  altitudes  the  maximum  temperature  exceeds  1800  F. 

The  boundaries  are  defined  by  laminar  heating  at  the  forward  sta¬ 
tions  during  the  initial  portion  of  entry.  The  break  in  the  ini¬ 
tial  slope  of  the  boundary  is  due  to  the  onset  and  realization  of 
turbulent  heating  at  the  aft  stations. 

An  inequality  constraint  is  also  imposed  on  the  lift  coeffi¬ 
cient  which  is  limited  to  CL  =  0.535. 
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VARIATIONAL  TREATMENT 


Unlike  the  early  energy  approximation  studies  where  only  one 
differential  equation  had  to  be  analyzed,  the  present  investigation 
involves  four  differential  equations  of  motion  and  two  inequality 
constraints  on  the  control  and  state  variables.  This  more  complex 
problem  requires  the  application  of  optimal  control  theory  (Ref.  26) 
to  formulate  the  necessary  conditions  for  a  maximum  crossrange  turn¬ 
ing  re-entry.  These  conditions,  together  with  the  state  equations 
(11)  to  (14),  form  a  two-point  boundary  value  problem  that  must  be 
solved  numerically/iteratively  for  a  specific  set  of  boundary  con¬ 
ditions  . 


Necessary  Conditions 

The  Hamiltonian  for  the  unconstrained  problem  is 


H 


DV  /^o^  s^n  ^  v  \ 

W  +  \  (  — WV - rSinx  tan  ®) 


X  (2  cos  x)  +  v 

\r  v  A  \r  cos  tp/ 


(20) 


and  the  augmented  Hamiltonian  for  the  constrained  problem  is 

8  '  H  +  VUBCUB  +  VMLCML  (21) 

where  A£,  A  ,  A^,  A^  are  undetermined  multipliers,  Cyg  m  h  -  hy^  >  0 

is  the  underbody  heating  constraint,  »  ^max  "  ^  >  0  is  the 

maximum  lift  constraint,  hyy  is  given  by  the  thermal  boundaries 
in  Fig.  3,  L  is  the  maximum  lift  with  maximum  CT  *  0.535, 

tn^x  L 

and  v^,  are  undetermined  constraint  multipliers  that  must 

satisfy  the  following  conditions 
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max  — 

The  two  additional  terms  in  the  augmented  Hamiltonian  correspond 
to  the  treatment  of  inequality  constraints  bv  the  technique  of 
Valentine  (Refs.  27  and  28). 

Substitution  of  Eqs.  (10),  (16),  (18),  (19)  into  Eqs.  (20), 

'Y/ 

(21)  and  differentiation  of  H  with  respect  to  the  four  state 
variables  leads  to  the  following  Euler  equations: 
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where,  by  differentiating  Eq.  (9), 


The  partial  derivatives  dCyg/dE  and  dC^/dE  are  given  in  Appen 
dix  C.  Since  C^g  and  are  not  functions  of  x>  CP>  or  A, 

their  partial  derivatives  do  not  appear  in  Eqs.  (25)  to  (27). 

For  optimality,  the  Hamiltonian  must  be  stationary  with  re¬ 
spect  to  the  two  control  variables.  The  Hamiltonian,  for  the  un¬ 
constrained  case,  or  unconstrained  portion  of  a  trajectory,  can  be 


written  as 


where 


2 

H  =  A  +  B  sec  m-  +  C  tan  p 


YD  ^  *  \  7  Ca"  Sln  * 


.  V  .  V  sin  x 

A  —  cos  x  +  “ 

<P  r  ^  A  r  cos  cp 


,  „  2WV  8oro 

ECD  ,  ~r 

C?  pSgo  Vr 
L 


goro  V 


Bank  angle  is  determined  by  the  partial  derivative  dH/dp  ■  0 
and  is  obtained  by  differentiating  Eq.  (29).  The  result  is 


i. 

L 

K 

L 

L 

V. 

0 

i . 

I 

n 

r. 

i: 

i. 


m 

i 


tan  m. 


_C_ 

2B 


Substitution  of  Eq.  (33)  into  Eq.  (25)  results  in 


H 


A  +  B  - 


(33) 


(34) 


which  is  now  a  function  of  only  one  variable,  r.  The  value  of  r 
which  minimizes  the  Hamiltonian  [Eq.  (34)  ]  is  determined  by  a  one 
dimensional  search. 


The  Hamiltonian  as  given  by  Eq.  (29)  assumes  a  parabolic  drag 
relationship;  however,  a  similar  investigation  can  be  carried  out 
for  a  fourth-order  drag  model.  The  Hamiltonian,  Euler  equations, 
and  optimum  bank  angle  for  this  more  realistic  drag  model  are  given 
in  Appendix  D. 

For  the  case  involving  the  underbody  heating  constraint,  the 
Hamiltonian  H  must  also  be  stationary  with  respect  to  the  control 
variables.  This  results  in 


SH  .  SH  ,  v  ^CUB  _  Q 

Sp  Sp  UB  dp. 


(35) 


SH  SH  ,  ^CUB 
&  "  +  VUB  ~^r" 


0 


(36) 


Similarly,  for  the  maximum  lift  constraint  portion  of  a  trajectory 


SH  _  SH  5CML  ,  0 

£p  Sp  ML  dp 


(37) 


o>H  =  oH  .  v  ^StL 
5r  Sr  ML  Sr 


0 


(38) 
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Eq.  (41)  was  used  as  a  check  on  the  accuracy  of  the  numerical  pro¬ 
cedure.  As  previously  mentioned,  the  initial  conditions  of  E,  y, 
<P,  A,  t  are  specified  as  well  as  the  final  value  of  E;  final 
values  of  y,  9,  A,  t  are  open.  The  integration  forward  in  time 
can  proceed  once  appropriate  choices  are  made  for  the  initial 
values  of  the  undetermined  multipliers.  One  of  these  multipliers 
is  selected  by  the  scaling  condition,  i.e.,  A  *  -  1;  a  second 
multiplier  can  be  determined  from  Eqs.  (29)  to  (32)  and  (41), 

which  at  t  reduces  to 
o 


4B  B 


(42) 


Since  A.  is  a  function  of  r,  the  Hamiltonian  har  to  be  mini- 
o 

mized  before  computing  A^  from  Eq.  (42).  A  family  of  trajec¬ 
tories  can  be  generated  in  terms  of  the  remaining  two  multiplier 

initial  values  A-  and  A  .  A  point  on  the  footprint  is  ob- 

o 

tained  by  selecting  A  and  iteratively  searching  for  that  value 

Ko 

of  A<p  which  results  in  a  solution  that  satisfies  all  necessary 

conditions  and  for  which  the  final  value  of  A  is  zero  (trans- 

f 

versality  condition) .  A  neighboring  point  on  the  footprint  is  ob¬ 
tained  by  varying  A  and  A  ,  utilizing  continuity  properties 

*o  ^  o 

of  the  state  and  multiplier  variables  with  respect  to  boundary  con¬ 
ditions. 


The  constraint  multipliers  for  the  constrained  portion  of  the 
trajectory  are  obtained  from  Eqs.  (36)  and  (38),  resulting  in 


UB 


ML 


dH 

dr 


dH 

d 


/ 


3cub 

(43) 

^°ML 

(44) 
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For  the  case  where  both  constraints  apply  at  the  same  time, 
Eqs.  (39)  and  (40)  are  solved  simultaneously 


d°ML 

-1 

dH 

dp 

dp 

dp 

SCUB 

d°ML 

dH 

dr 

dr 

Sr. 

(45) 


At  each  time  interval  a  one  dimensional  search  is  made  for 
the  altitude  that  minimizes  the  Hamiltonian.  For  each  altitude 
during  the  search  Eqs.  (33),  (10),  (16),  and  (17)  are  used  to  com¬ 
pute  the  unconstrained  optimum  bank  angle  and  angle  of  attack,  and 
the  results  compared  with  the  heating  constraint  given  by  Fig.  3. 
If  the  heating  constraint  is  violated  the  unconstrained  numerical 
results  are  discarded  in  favor  of  the  nearest  angle  of  attack  that 
will  satisfy  the  constraint.  This  is  accomplished  by  the  follow¬ 
ing  iterative  procedure,  initiated  with  the  numerical  value  of  the 
unconstrained  angle  of  attack: 


2. 


a, 


a.  ,  +  Act. 
l-l  l 


3-  hUB1  "  him(V  v> 

4.  test  |h  -  h  |  <  « 


where  c  is  an  arbitrary  small  number,  hTTD  is  a  double  table 
look-up  with  data  taken  from  Fig.  3,  and  dh^/da  is  a  piecewise- 
constant  function  taken  from  the  data  stored  in  the  double  table 
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look-up.  The  four  steps  are  repeated  until  the  altitude  difference 
is  less  than  €.  The  cosine  of  the  bank  angle  is  then  computed 
from  Eqs.  (10)  and  (15)  to  (17),  i.e.. 


cos  u 


2 

WV  fo  /foN  V 
g  L  V  \rJ  ~  r 


(46) 


If  cos  u  >  1,  an  impossible  equilibrium  flight  condition  exists 
for  that  altitude  and  the  next  value  of  attitude  is  selected  in 
the  search  to  minimize  H.  If  cos  p.  >  1  for  all  values  of  alti¬ 
tude,  the  heating  constraint  cannot  be  satisfied  and  the  run  is 
terminated. 

The  maximum  lift  coefficient  constraint  is  satisfied  by  com¬ 
puting  lift  from  Eq.  (10),  maximum  lift  from  Eq.  (16)  using 

CT  *  CT  ,  and  verifying  that  L  <  L  .  If  L  >  L  ,  the  un- 
l  l  max  max 

max 

constrained  values  of  bank  angle  and  angle  of  attack  are  discarded 
and,  instead,  Eq.  (46)  is  used  to  compute  the  constrained  bank 
angle  with  L  «  L  ,  and,  the  search  is  continued  for  the  value 
of  r  that  minimizes  H. 
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NUMERICAL  RESULTS 


The  numerical  data  employed  for  the  re-entry  vehicle  corre¬ 
sponds  to  one  of  the  early  (1970)  space  shuttle  orbiter  configura¬ 
tions.  A  FORTRAN  IV  digital  computer  program  was  used  to  calculate 
iteratively  bo^h  unconstrained  and  constrained  optimum  trajector¬ 
ies.  The  footprints  for  both  are  shown  in  Fig.  4  as  a  locus  of 
terminal  latitudinal  and  longitudinal  points  for  a  fixed  terminal 
energy  state.  The  unconstrained  footprint  was  traced  from  the 
easily  determined  maximum  downrange  point  (<p  *=  X  *  M-  “  *  \  *  °)> 

T  /\ 

back  along  the  north  branch  and  past  the  maximum  crossrange  point. 

No  attempt  was  made  to  define  the  rear  of  the  footprint,  which  is 
believed  to  involve  zigzagging  paths.  The  effect  of  the  heating 
constraint  is  to  shrink  the  footprint,  as  would  be  expected,  with 
the  maximum  downrange  and  crossrange  reduced  by  64  and  450 
nautical  miles,  respectively. 

Figures  5,  6,  and  7  show  the  trajectory  characteristics  for 
typical  high  and  low  crossrange  solutions,  with  and  without  heating 
constraints.  For  the  sequence  of  runs  that  include  the  heating 
constraint,  the  trajectories  are  initially  constrained,  remain  so 
most  of  the  time,  and  thereafter  continue  to  be  unconstrained. 


The  high  sensitivity  of  the  final  state  conditions  with  re¬ 
spect  to  the  initial  conditions  of  the  multipliers  is  demonstrated 
in  Fig.  8.  Optimal  runs  have  zero  bank  at  the  terminal  conditions. 
This  can  be  shown  by  substituting  Eqs.  (31)  and  (32)  into  Eq.  (33) 


tan  p. 


^  psg2 

Y  &Q 

r  _  2 


4Xewcd  , 

goro  V2 

2  "  r 

c? 

(47) 
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Fig.  5  Typical  Reduced-Order  Trajectories  for  Space  Shuttle 
Re-entry  -  Latitude  vs  Longitude 
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Fig.  6  Typical  Reduced-Order  Solutions  for  Space  Shuttle  Re-entry 
Altitude  vs  Velocity 
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and  noting  that  the  transversality  condition  calls  for  \  ■  0. 

Neighboring  solutions  have  hard  right  or  hard  left  maneuvers  de¬ 
pending  upon  the  sign  of  ^  at  the  specified  value  of 

■  -  20,800,000  feet.  Fortunately,  in  spite  of  this  sensi¬ 
tivity,  it  is  clear-cut  how  to  change  iteratively  the  numerical 
value  of  to  achieve  =0. 

The  maximum  lift  constraint  was  added  to  the  computer  program 
in  the  anticipation  that  large  bank  angles  at  high  altitudes,  or 
high  altitudes  resulting  from  the  heating  constraint,  would  call 
for  excessive  lift  capability.  However,  this  never  occurred.  In¬ 
stead,  it  was  noted  that  the  lift  constraint  was  necessary  at  the 
end  of  some  of  those  trajectories  that  involved  hard  turning  maneu¬ 
vers.  For  those  trajectories  with  -  0,  i.e.,  Uf  *  0,  the 

lift  constraint  was  never  required. 

In  Appendix  B  it  is  shown  analytically  that  the  optimum  angle 
of  attack  for  the  unconstrained  case,  or  unconstrained  portion  of 
a  trajectory,  is  that  which  maximizes  L/D.  Substitution  of  the 
numerical  data  into  Eqs.  (B-22)  and  (17)  produces  an  angle  of 
attack  of  15.6°.  No  attempt  was  made  to  capitalize  on  this  ana¬ 
lytical  result  when  writing  the  computer  program;  however,  all  of 
the  numerical  solutions  for  the  unconstrained  cases,  as  well  as  the 
unconstrained  portions  of  the  constrained  cases  (see  Fig.  9),  con¬ 
sist  of  angle  of  attack  time  histories  that  remains  constant  at 
15.6  ,  thus  verifying  the  analytical  finding.  It  has  been  shown 
in  the  classical  literature  that  for  two  dimensional  flight  over  a 
flat  earth  maximum  L/D  will  provide  the  most  range.  The  analytic 
solution  of  Appendix  B  and  all  of  the  numerical  results  of  this  re¬ 
port  (without  constraints)  show  that  maximum  L/D  is  also  the  op¬ 
timum  control  that  maximizes  range  and  crossrange  for  three  dimen¬ 
sional  turning  flight  over  a  spherical  earth  with  an  inverse -square 
gravity  field. 
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Fig.  9  Typical  Angle  of  Attack  Time  Histories  for  Reduced-Order 
Solutions 


Throughout  the  investigation  and  computer  programming  no  de¬ 
liberate  effort  was  made  to  minimize  machine  time.  However,  there 
are  several  worthwhile  approaches  that  could  be  exploited.  The 
first  is  to  take  advantage  of  the  analytic  solution,  maximum  L/D 
angle  of  attack,  for  the  unconstrained  re-entry  and  unconstrained 
portions  of  the  trajectories.  The  optimal  controls  are  given  ex¬ 
plicitly  in  Appendix  B  by  Eqs .  (B-22)  and  (B-24)  .  No  numerical 
search  is  necessary  to  minimize  the  Hamiltonian;  however,  a  simple 
iterative  procedure  may  be  required  to  determine  the  altitude  from 
Eq.  (B-7) .  The  second  approach  at  streamlining  the  computer  pro¬ 
gram  would  be  to  utilize  the  constants  of  motion  from  Noether's 
theorem  (Appendix  E) ,  together  with  the  constant,  zero  value  of 
the  Hamiltonian.  These  four  constants,  given  by  Eqs.  (C-l) ,  (C-2) , 
(C-3) ,  and  (41),  eliminate  the  need  to  numerically  integrate 
Eqs.  (24)  to  (26).  A  third  approach  would  be  to  utilize  a  system 
of  differential  equations  where  energy  is  used  as  the  independent 
variable  (Appendix  B) .  This  would  eliminate  the  need  to  numeri¬ 
cally  integrate  Eq.  (11).  Although  the  Hamiltonian  would  no  longer 
be  constant,  this  would  be  compensated  by  the  elimination  of  the 
energy  multiplier,  ^ . 
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Ij  APPENDIX  A 

DERIVATION  OF  EQUATIONS  OF  MOTION 

The  flight  path  equations  of  motion  could  be  derived  from 
first  principles;  however,  a  rigorous  development  would  be  rather 
lengthy  and  not  in  keeping  with  the  purpose  of  this  report.  In¬ 
stead,  utilization  will  be  made  of  Ref.  29  which  provides  a  thorough 
derivation  and  requires  only  a  few  additional  substitutions  to  ob¬ 
tain  Eqs.  (1)  to  (6).  Since  many  of  the  variables  of  Ref.  29  have 
definitions  and  symbols  different  from  those  of  this  report,  only 
Appendix  A  will  employ  the  symbols  of  Ref.  29. 

The  curvilinear  coordinates  X  and  Y  are  related  to  longi¬ 
tude  t  and  latitude  ^  (Fig.  A-l)  by  simple  arc  length  formulae 
given  by  Eqs.  (6)  on  p.  61  of  Ref.  29 


X  =  r  r  ,  Y  =  r  A  (A-l) 

o  o 

where  the  constant  r  is  the  radius  of  the  earth.  The  resulting 

o 

curvilinear  velocities  are 


| 

i 


•  •  • 

X  =  r  t  ,  Y  =  r  *  (A-2) 

o  o 

The  angular  velocities  of  the  wind  axes  (p  ,  q.  r  )  centered  at 

W  W  W 

the  vehicle's  eg  (Fig.  A-2)  are  related  to  earth  fixed  axes  cen¬ 
tered  at  some  point  on  the  surface  of  the  earth  (Fig.  A-l)  by 
Eq.  (17)  on  p.  63.  Substitution  of  (A-l)  and  (A-2)  into  this  equa 
tion  results  in 


Pw 

1  . 

s* 

cs 

- 

r 

<•  m 

w 

1  0  -  sin  y 

0  cos  \x  sin  m.  cos  y 

0  -  sin  |i  cos  p  cos  y 


* 


28 


Vehicle  x 


Vehicle’s 

Parallel 


Equator 


Vehicle's 

Meridian 


Fig.  A-l  Coordinate  Systems  for  Flight  Over  Spherical  Earth  (Ref.  29) 


Fig.  A-2  System  of  Rotations 


from  Local  Horizon  to  Wind  Axes  (Ref.  29) 


cos  y  cos  x 

sin  m-  sin  y  cos  x 
-  cos  |i  sin  x 

cos  \i  sin  y  cos  x 
+  sin  m-  sin  x 


cos  y  sin  x 

sin  Li  sin  y  sin  x 
+  cos  |i  cos  x 

cos  ^  sin  y  sin  x 
-  sin  p.  cos  x 


-  sin  7 
sin  p.  cos  7 

cos  Li  cos  7 


X 

T  COS  x 

t  sin  X 


(A-3) 


where  x  is  the  azimuth  heading  angle  (0°  due  west,  90°  due 
north,  etc.),  7  is  the  flight  path  angle  above  horizontal,  and 
U  is  the  bank  angle  (positive  right  wing  down) .  The  kinematic  re¬ 
lationships,  which  when  integrated  defines  the  vehicle's  position 
with  respect  to  the  earth,  are  given  by  Eqs.  (22)  on  p.  64.  Sub¬ 
stitution  of  (A-l)  and  (A-2)  into  this  equation  results  in 


V  cos  y  cos  X 
r  cos  X 


(A-4) 


;  v 

X  ■  —  cos 
r 


7  sin  x 


(A-5) 


h  =  V  sin  7 


(A-6) 


where  h  is  the  altitude  and  r  -  r  +  h  is  the  distance  of  the 

o 

vehicle  from  the  center  of  the  earth.  Gravity,  the  usual  inverse- 
square  law,  is  given  by  Eq.  (25)  on  p.  65 

2 

8  ■  go  (r"+h)  (w) 

o 

The  dynamical  relationships  are  given  by  Eqs.  (26)  on  p.  65.  Sub¬ 
stitution  of  (A-7)  into  this  equation  and  assuming  no  thrust,  a 
coordinated  turn  with  no  side  force,  and  a  nonrotating  earth  re¬ 
sults  in 
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V  ■  -  (D/m)  -  g  sin  y 

CA-8) 

rw  "  (g/V)sin  p  cos  y 

(A-9) 

*  (L/mV)  -  (g/V)cos  p  cos  y 

(A-10) 

Substitution  of  rw  and  qw  from  (A-3)  into  (A-9)  and 
becomes 

(A-10) 

•  . 

cos  p  cos  y  -  y  sin  p  +  (V/r)cos  y  sin  p 

2 

+  (V/r)cos  y  cos  y  cos  p  tan  A  «  (g/V)sin 

(A-ll) 

u  cos  y 

• 

X  sin  p  cos  y  +  y  cos  p  -  (V/r)cos  y  cos  pi 

2 

+  (V/r)cos  y  cos  y  sin  p  tan  A 

(A-12) 

■  (L/mV)  -  (g/V) cos  p  cos  y 

Multiplying  (A-ll)  by  cos  u  and  (A-12)  by  sin  u  and  then  adding, 
and  similarly  multiplying  (A-ll)  by  -  sin  u  and  (A-12)  by  cos  u 
and  then  adding,  results  in 

• 

X  cos  y  +  (V/r)cos2y  cos  y  tan  A  «  (L/mV)sin  p 

(A-13) 

• 

y  -  (V/r)cos  y  =  (L/mV)cos  p  -  (g/V)cos  y 

(A- 14) 

The  six  degree  of  freedom  equations  are  composed  of  (A-4) ,  (A-5) , 

(A-6),  (A-8),  (A-13) ,  and  (A- 14) .  Since  h  =  r,  these  equations 
can  be  summarized  as; 
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(A— 15) 


V 

7 


x 

r 


T 


"  m  “  8  sln  7 


2v  cos 


g  ,  V 

-  cos  y  +  —  cos  7 


L  sin  (j. 
mV  cos  7 


-  ^  cos  7  cos  \  tan  X 


V  sin  7 


V  cos  7  cos  Y 
r  cos  X 


(A- 16) 

(A-17) 
(A- 18) 
(A-19) 


*  V 

^  ■  —  cos  7  sin  v  (A- 20) 

The  differences  between  the  above  six  equations  and  Eqs.  (1)  to 
(6)  are  due  to  differences  in  definitions  and  symbols.  The  longi¬ 
tude  t  increases  heading  west  (Fig.  A-l) ,  whereas  A  in  (6) 
increases  heading  east;  the  bank  angle  p  above  is  positive  right 
wing  down,  whereas  in  (2)  and  (3)  it  is  positive  left  wing  down; 
and  the  azimuth  heading  angle  in  (3),  (5),  and  (6)  is  0°  due 
south,  90°  due  east,  etc.  If  the  following  four  changes  are 
made,  Eqs.  (1)  to  (6)  will  be  obtained 

T  — *  “A  ,  IL  -»  “P. 

*  <P  ,  X  ->  270°  -  y 


An  almost  identical  system  of  equations  of  motion  can  also  be 
found  in  Ref.  30. 
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APPENDIX  B 


CHANGE  OF  INDEPENDENT  VARIABLE  FROM  TIME  TO  ENERGY 


Based  on  the  somewhat  limited  numerical  results  of  this  report 
it  is  possible  to  convert  the  independent  variable  from  time  to 
velocity,  or  heading  angle,  or  latitude,  or  longitude,  since  all  of 
these  state  variables  have  been  noted  to  vary  monotonically  for  the 
reduced-order  solutions.  However,  energy  is  the  only  state  variable 
which  is  certain  always  to  vary  monotonically  for  all  unpowered 
flight  problems,  even  for  the  complete  six  degree  of  freedom  equa¬ 
tions.  From  Eqs.  (11)  and  (1)  to  (6) 


dt  _W_ 
dE  “  “  DV 


(B-l) 


dy  dy  dt 
dE  "  dt  dE 


Lgo  ,  Wg  W 

— s'  cos  M-  +  — J  cos  y  -  r—  cos  7 

DV  DV  r 


(B-2) 


dy  ^  dy  dt 

dE  "  dt  dE 


sin  ii  u 

- +  r—  cos  7  sin  y  tan  <p 

cos  7  Dr  ' 


(B-3 


dr  _  dr  dt_  W 
dE  "  dt  dE  "  “  D 


(B-4) 


d«P  dtp  dt  W 

dE  “  dF  dE  =  D?  C°S  7  C°S  X 


(B-5) 


dA  dA  dt  _  _W_  cos  7  sin  y 

dE  dt  dE  Dr  cos  <p  '  ”  ' 

Utilizing  the  same  previous  assumptions,  that  the  flight  path 
angle  7  is  small  and  the  angular  rate  7  is  zero,  (B-2)  and  (B-4) 
are  converted  to  algebraic  expressions 
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COS  p 


(B-7) 


W  L  8or 


D  D  ,  ,.2v 

(gr  -  V  ) 


r  =  r  +  h 
o 


(B-8) 


The  system  of  differential  equations  is  reduced  to  third  order: 


dv  l 
dE  =  D 


g0  go 

- r—  cos  p  sin  y  tan  9  -  —zr  sin  p 

(gr  -  V2)  V 


(B“9) 


L  8o 
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dA  =  L  °o  cos  p  sin  y 
dE  ”  "  D 


(gr  -  V2)  cos  ” 


(B-ll) 


The  coefficients  involving  gravity  and  velocity  are  converted  with 
the  use  of  the  energy  equation  (9) 
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Equations  (B-9)  to  (B-ll)  become 
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dtp  L  cos  4 
dE  "  "  D  r2 

2E  +  -f 

dA  l  cos  4  sin  x 
dE  "  D  r2  cos  <p 

2E  +  ~ 
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and  the  Hamiltonian  can  be  written: 


H  =  * 


cos  4  /sin  \ 

^2  \cos  cp  )  A 


-  ^  sin  cpj.  -  cos  y) 


2E  +  -2- 
r 


(B-17) 


sin  4 


r2  x 

2(E  +  f) 


Although  there  are  three  control  variables  in  (B-17),  namely  r,  4, 
and  a  (L/D  is  a  function  only  of  a) ,  r  plays  a  very  insignifi¬ 
cant  part,  almost  like  a  constant.  This  is  vastly  different  from 
the  Hamiltonian  given  by  Eqs.  (29)  to  (32),  where  the  air  density  p 
which  is  strongly  dependent  on  r,  appears  in  the  A  and  B  terms 
For  example,  at  a  typical  initial  condition  of  V  *=  24,100  ft/sec 
and  h  *  233,000  ft,  a  rather  large  change  in  r  of  10,000  ft 
results  in  a  52%  change  in  p  and  only  a  0.38%  change  in 
(2E  +  rQ/r)  and  a  0.11%  change  in  2(E  +  rQ/r)  .  Advantage  can 
be  taken  of  this  property  by  employing  the  following  different  pro- 
cedure  for  minimizing  the  Hamiltonian  while  satisfying  the  y  ■  0 
condition,  i.e.,  (B-7) . 


Equation  (B-17)  can  also  be  written  as 


H  ■  -*[A  cos  p  +  B  sin  p] 


where 


r2  NCOS  <p  )  A 


X  sin  <p 
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*  \  cos  x) 
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(B-18) 


(B-19) 


(B-20) 


First,  a  reasonable  guess  for  r  would  be  made;  e.g.,  the  value 
of  r  from  the  previous  integration  interval.  Then  the  Hamiltonian 
is  minimized  with  respect  to  a,  which  results  in  maximizing  L/D. 
For  drag  polars,  as  given  by  Eq.  (18),  it  can  be  shown  that 
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The  Hamiltonian  is  then  minimized  with  respect  to  p 


dH 

dp 


It[- 

Dl 


A  sin  p  +  B  cos  p]  =  0 


Consequently 


tan  p 


B 

A 


(B-23) 


(B-24) 


The  7*0  condition  is  then  satisfied  by  iteratively  computing  a 
new  value  of  r  from  (B-7)  .  This  procedure  of  computing  A,  B,  ijl, 
and  r  can  be  repeated,  if  necessary,  until  the  change  of  r  is 
less  than  some  small  specified  value. 


APPENDIX  C 


DERIVATION  OF  PARTIAL  DERIVATIVES 

The  partial  derivatives  of  the  augmented  Hamiltonian  H  with 
respect  to  the  four  state  variables  E,  y,  <P»  A  are  easily  derived 
and  result  in  the  four  Euler  equations,  (24)  to  (27).  All  of  the 
other  derivatives  required  for  the  computer  program  are  obtained  by 
partial  differentiation  of  H,  with  respect  to  E,  R,  p. 

Rules  of  partial  differentiation,  concerning  which  quantites  are 
fixed  or  variable,  are  the  same  as  those  governing  the  development 
of  the  Euler  equations. 

S(CUB’  W/5E 

Partial  derivatives  taken  with  respect  to  E  are  carried  out 
for  fixed  values  of  the  control  variables  r,  p,  and  the  state 
variables  y,  cp,  A.  All  other  variables,  such  as  V,  a,  C L,  L,  etc. 
are  not  fixed.  Repeating  Eq.  (22), 

CUB  -  h  -  hUB  <«*» 

where  h  =  r  -  rQ,  and  h^  is  given  by  the  thermal  boundaries 
in  Fig.  3.  Differentiating  (C-l)  with  respect  to  E, 

^CUB  =  ^hUB  BV  ^UB  Ba_  BV  (r 

Be-  bv  be  B:*  bcl  bv  Se  k  } 

where  BhITT1/BV  and  Bh.m/Ba  are  piecewise  constant  functions  com- 

Ud  Ud 

puted  from  h._  data  stored  in  the  double  table  look-up.  From  the 
Ud 

energy  equation  (9), 


From  the  angle  of  attack  equation  (17), 


5^-  -98  sec2(1.33  CL) 
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Combining  Eqs .  (10)  and  (16),  and  differentiating. 
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Repeating  Eq.  (23), 
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where  l>max  and  L  are  obtained  from  Eqs.  (10)  and  (16) 
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Differentiating  (C-7)  to  (C-9)  with  respect  to  E, 
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Partial  derivatives  taken  with  respect  to  r  are  carried  out 
for  fixed  values  of  the  multipliers,  the  state  variables  E,  x» 

A,  and  the  control  variable  p.  All  other  variables,  such  as  V, 
a,  C^,  L,  etc.,  are  not  fixed.  The  Hamiltonian  is  given  by 
Eq.  (29) 

H  =  A  +  B  sec2p  +  C  tan  p  (C-13) 


where  A,  B,  C  are  given  by  Eqs.  (30)  to  (32).  Differentiating 
(C-13)  with  respect  to  r, 
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2 
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+  fftan 
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Since  A,  B,  C,  as  well  as  CT  ,  L  ,  L  that  follow,  are  functions 

L  max 

of  r  and  V,  the  partial  derivatives  are  derived  by 


Sa  /  San  .  (  sa\  Sv 

Sr  "  \Sr/ ,  +  \bv)  Sr 

V  r 


(C-15) 


where  the  subscripts  indicate  which  variables  are  fixed.  From 
Eqs.  (9)  and  (15) 
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Thus, 
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The  derivative  SClro/Sr  is  obtained  by  differentiating  Eq.  (Ol) 

Ud 


SCUB  -  .  ^UB  SV  .  ahUB  Sa  5CL 

Sr  SV  or  Sa  SCL  Sr 


(021) 


where,  by  differentiating  Eq.  (05) 
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All  other  derivatives  in  (021)  have  been  previously  defined. 
The  derivative  SC^/Sr  is  obtained  from  Eq.  (07) 
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Sr  Sr 

where,  by  differentiating  Eqs .  (08)  and  (09) 
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(C-25) 


dL 

dr 


wv2 


2 

gQr  cos  p 


C^,  C^/dp 

Partial  derivatives  taken  with  respect  to  p  are  carried  out 
for  fixed  values  of  the  multipliers,  the  state  variables  E,  x> 

A,  and  the  control  variable  r.  Since  V  is  a  function  only  of 
E  and  r,  both  of  which  are  fixed,  V  also  is  fixed.  All  other 
variables,  such  as  a,  C^,  L,  etc.,  are  not  fixed.  Differentiating 
Eq.  (013)  with  respect  to  p, 

=  [2B  tan  p  +  C]sec2p  (026) 


Differentiating  Eq.  (Ol)  with  respect  to  p 
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UB  da 
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da  dC^  dp 


(027) 


where  the  first  two  derivatives  have  been  previously  defined  and, 
by  differentiating  Eq.  (05) 
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Combining  Eqs .  (C-7)  to  (09)  and  differentiating  with  respect 
to  u  results  in 
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APPENDIX  D 


DRAG  POLARS  OF  FOURTH  ORDER 

The  mathematical  formulation  and  computer  program  described  in 
this  report  applies  only  to  those  cases  where  the  aerodynamic  data 
are  represented  by  parabolic  drag  relationships,  as  given  by 
Eq.  (18).  Because  of  the  uncertainty  in  the  aerodynamic  data,  par¬ 
ticularly  during  the  preliminary  design  phase,  this  simple  drag 
polar  should  be  adequate  for  many  studies.  However,  as  more  wind 
tunnel  data  become  available,  it  may  be  necessary  to  improve  the 
polar  curve  fit  by  adding  a  fourth-order  term 


CD*CD  +CD,CL+CD  X 

0  4 


(D-l) 


where  C_  is  always  positive,  both  for  straight  and  swept  wing 

c4 

re-entry  vehicles.  In  the  following  analysis,  this  is  not  only 
important,  but  fortunate.  For  this  fourth-order  drag  polar,  the 
Hamiltonian  is 


2  4 

H  =  A  +  B  sec  ^  +  C  tan  p  +  F  sec  p 


(D-2) 


where  A,  B,  and  C  are  given  by  Eqs .  (30),  (31),  (32),  and 
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The  Euler  equations  for  A  ,  *  ,  A.,  as  given  by  Eqs.  (25)  to  (27) 

)(  T 

do  not  change;  however 
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where,  as  before. 


Bank  angle  is  determined  from  SH/Sm.  =  0  and  results  in  the 
following  cubic  equation: 
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Thus  tan  4  +  a  tan  p.  +  b  ■  0.  The  roots  of  this  cubic  equation 

depend  on  the  sign  of  (b2/4  +  a^/27)  .  It  can  be  shown  that  this 

quantity  is  always  positive,  providing  both  Cn  and  Cn  have 

D  2  D  4 

CL  CL 

the  same  sign,  and  therefore  there  is  only  one  real  root  (the  other 
two  are  imaginary  and  can  be  ignored)  .  Since  both  of  these  coeffi¬ 
cients  are  positive,  there  is  only  one  solution,  namely 


where  a  and  b  are  given  by  Eqs .  (D-6)  and  (D-7) .  Substitution 
of  Eq.  (D-8)  into  (D-2)  results  in 

H  =  A  +  B  +  CG+F  +  BG2+  2FG2  +  FG4 


which,  like  Eq.  (34),  is  a  function  of  only  one  control  variable, 
namely,  r. 


(D-9) 


APPENDIX  E 


CONSTANTS  OF  THE  MOTION  FROM  NOETHER'S  THEOREM 


Noether's  theorem  (Ref.  31)  applies  to  variational  problems 
that  possess  symmetry.  Symmetry  exists  if  one  can  subject  a  system 
of  equations  of  motion  to  a  certain  operation  (transformation)  and 
the  system  appears  exactly  the  same  after  the  operation.  For  ex¬ 
ample,  Noether's  theorem  shows  that  if  the  coordinate  system  that 
defines  the  state  variables  can  be  translated,  rotated,  or  shifted 
in  time  without  changing  the  dynamics  of  the  problem  then  for  each 
transformation  a  constant  of  the  motion  can  be  derived.  This  ap¬ 
plies  not  only  to  the  variational  formulations  of  classical  and 
quantum  mechanics,  for  which  the  theorem  was  originally  derived 
(Ref.  31),  but  also  to  modem  optimal  control  problems  (Ref.  32), 
for  which  the  constants  of  the  motion  involve  both  state  and  multi¬ 
plier  variables. 

The  coordinate  system  for  Eqs.  (1)  to  (6)  and  the  reduced  sys¬ 
tem  (11)  to  (14)  are  such  that  rotations  about  the  three  axes  of 
the  reference  system  which  define  A,  <p,  and  \  can  be  made  with¬ 
out  affecting  either  the  dynamics  or  the  optimal  solutions.  This 
would  not  be  true  for  the  case  with  earth  rotation,  since  a  change 
in  orientation  of  the  earth's  axis  of  rotation  would  result  in 
different  optimal  trajectory  solutions. 

Using  Noether's  theorem,  H.  G.  Moyer  of  the  Grumman  Research 
Department  derived  the  following  three  constants: 

C,  =  A  cos  A  sec  cp  -  X  cos  A  tan  qp  +  X  sin  A  (E-l) 

lX  A  qp 

C0  ■  *  sin  A  sec  qp  -  X  sin  A  tan  cp  -  X  cos  A  (E-2) 

2  x  4  cp  ' 

C3  -  \  <E-3) 
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Writing  these  three  equations  in  matrix  form  and  inverting  the 
square  matrix,  the  multipliers  are  expressed  in  terms  of  latitude, 
longitude,  and  the  three  constants 


cos  A  cos  <p  +  C2  sin  A  cos  9  +  Cg  sin  9 


XA  "  C3 


A  «  C,  sin  A  -  C0  cos  A 
9  1  2 


(E-4) 


(E-5) 


(E-6) 


where  the  determinant  of  the  matrix  given  by  (E-l)  to  (E-3)  is 
sec  9,  which  is  singular  only  when  cp  -  ±90°  (crossrange  exceeds 
5000  nautical  miles)  . 

Equations  (E-l)  to  (E-6)  have  been  verified  analytically  for 
the  complete  system  (1)  to  (6) ,  the  reduced  system  (11)  to  (14) , 
and  the  reduced  system  (B-14)  to  (B-16)  with  energy  as  the  inde¬ 
pendent  variable. 

The  significance  of  the  three  constants  of  the  motion  is  two¬ 
fold.  First,  it  eliminates  the  step-by-step  numerical  integration 
of  A^  and  A^.  Since  A^  =  C^,  Eqs.  (B-17) ,  (B-19),  (B-20)  and 
(B-23)  can  be  expressed  in  terms  of  the  constants  and  state  vari¬ 
ables,  thereby  completely  eliminating  the  use  of  multipliers  for 
the  energy-independent-variable  system.  Similarly,  if  advantage 
is  taken  of  the  fact  that  the  Hamiltonian  is  constantly  zero,  the 
use  of  the  multipliers  for  the  time -independent -variable  system  is 
also  eliminated.  The  difficulty  in  determining  iteratively  the 
numerical  values  of  C^,  C2,  is  no  greater  than  determining 
the  initial  conditions  of  the  multipliers,  since 


■  — <  ,  I 


Care  must  be  taken,  however,  that  the  right  side  of  (E-4)  goes  to 
zero  at  the  end  of  each  optimal  solution.  Secondly,  since  the 
three  constants  are  valid  for  the  sixth-,  fourth-,  and  third-order 
systems,  as  well  as  problems  with  any  constraints  imposed,  it  is 
possible  to  predetermine  the  values  of  the  constants  as  a  function 
of  target  conditions  and  store  them  in  the  vehicle's  computer  to 
be  used  as  a  basis  for  an  on-board  guidance  scheme.  The  diffi¬ 
culty  here  is  that  the  constants  apply  to  a  nonrotating  earth  and 
may  have  to  be  modified  for  the  real  world  situation. 
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